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Abstract. We describe a new approach to heavy-quark fragmentation which is based on a resummed 
perturbative calculation and parametrization of power corrections, concentrating on the z — > 1 limit, 
where the heavy meson carries a large fraction of the momentum of the initial quark. It is shown that the 
leading power corrections in this region are controlled by the scale m(l — z). Renormalon analysis is then 
used to extend the perturbative treatment of soft and collinear radiation to the non-perturbative regime. 
Theoretical predictions are confronted with data on B-meson production in e + e~ annihilation. 



X 



PACS. 13. 66. Be Hadron production in e e 
- 12. 39. St Factorization 



interactions - 12.38.Cy Summation of perturbation theory 



1 Introduction 

The heavy-quark fragmentation function D(z,m 2 , /x 2 ) is 
the probability distribution to produce a heavy meson of 
a heavy quark. It depends on z, the momentum fraction 
of the meson, on the quark mass m 2 and on the factor- 
ization scale /j, 2 . The fragmentation function has a formal 
definition 2 as the Fourier transform 



D(z;ti 2 ) = 



1 



2nz 



— exp(ipy/z)F(py;n 2 ), (1) 

y- 



of the hadronic matrix element of a non-local operator on 
the light-cone (y 2 = 0): 



(2) 



F( P y;v 2 ) = 

' Tc{{0\#P{v)\H(p) + X)(H(p) + X\W(0)\0)^} 

x 



4 N r 



Here the final state is composed of the measured heavy 
meson (H) carrying momentum p plus anything else (X). 

We concentrate here on inclusive observables, the prime 
example being the single B-meson inclusive cross section 
in e + e~ annihilation, shown in Fig. Q] Cross sections of 
this sort can be written as a convolution 



da(x, Q 2 
dx 



dz 



C(x/z,Q 2 ]t , 2 )D(z; f , 2 ) (3) 



between a process-specific coefficient function C, describ- 
ing the hard interaction where the heavy quark is pro- 
duced, and the process-independent fragmentation func- 
tion D, defined in |T|I. describing the hadronization stage. 

* Invited talk (E.G.) at the HEP2003 Europhysics Confer- 
ence in Aachen, Germany. These proceedings are based on [Q. 




Fig. 1. ALEPH data for the energy distribution of produced B 
mesons at LEP1 compared to the next -to-leading order (NLO) 
calculation (full line). 

The common practice is to bridge the gap between ex- 
perimental data and fixed-order calculations in QCD by 
means of a fragmentation model, i.e. a given functional 
form for D(z; /i 2 ) in J2J with one or more free parameters. 
For heavy-quark fragmentation the most famous examples 
are [HE]- Upon excluding the more difficult z — > 1 re- 
gion such fits can indeed be performed. However, the gain 
is limited: the relation between the parameters in these 
models and the matrix elements cannot be made precise. 
The models provides no information about the underly- 
ing hadronization dynamics. Moreover, the universality of 
the extracted parameters is unclear. At large z fragmen- 
tation models simply fail to bridge the gap between the 
resummed perturbative calculation and the data. This has 
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been recently demonstrated in a clear way [H] by directly 
extracting the "non perturbative fragmentation compo- 
nent" from e + e~ data in moment space and then compar- 
ing the resulting distribution in z space to models. 

An important application of the heavy-quark fragmen- 
tation function which demonstrates these problems is in 
the description of B production in hadron colliders. The 
CDF collaboration found jo] an alarming discrepancy (a 
factor of 3) between the transverse-momentum distribu- 
tion of B + hadroproduction data and the standard treat- 
ment of this cross section, where a NLO calculation is 
convoluted with a Peterson model 0] for the fragmenta- 
tion function. In the latter the free parameter was set to a 
standard value based on e + e~ annihilation data. Ref. jZj 
applied a resummed perturbative calculation for the coef- 
ficient function and combined it with the relevant frag- 
mentation effect extracted from e + e~ data in moment 
space, concluding that the discrepancy is much smaller. 
This shows that the separation between the perturbative 
and non-perturbative ingredients of is very delicate. 
A naive application of (jSJ simply fails: if the perturbative 
ingredient C in © is taken at fixed order in a s , the re- 
quired "non-perturbative" ingredient D appears not to be 
the same in different processes. 

As heavy-quark production in hadron colliders becomes 
increasingly important experimentally, it is evermore ur- 
gent to correctly apply perturbative QCD to such cross 
sections, to separate in a systematic way between the 
perturbative and the non-perturbative ingredients, and fi- 
nally, to understand hadronization in a quantitative way. 
In particular, the parametrization of the fragmentation 
function D must eventually be understood in terms of its 
field theoretic definition QJ. 

Our approach to heavy-quark fragmentation is pri- 
marily a perturbative one: we start off with a perturba- 
tive calculation of the matrix element in J5J), replacing the 
outgoing meson by an on on-shell heavy quark, and treat 
non-perturbative effects, which make for the difference be- 
tween the quark and the meson, as corrections. Hadroniza- 
tion corrections are power-suppressed: they are inversely 
proportional to the mass of the heavy quark m. The per- 
turbative approach is appropriate so long as m 3> A. Thus 
it is definitely applicable to bottom, and probably, with 
some care, also to charm. 

It should be kept in mind that a perturbative calcula- 
tion is at all possible owing to two properties: (1) the pres- 
ence of the quark mass regulating collinear divergences; 
and (2) the inclusive nature of the observable, which guar- 
antees the cancellation of infrared singularities between 
real and virtual diagrams at any order in perturbation 
theory. This cancellation does leave, however, a signifi- 
cant trace in the expansion: Sudakov logarithms of 1 — z. 
This is why the 0(a s ) result shown in Fig. [T] diverges at 
z — ► 1, whereas the physical cross section vanishes at this 
limit. It is only upon summing the z — ► 1 singular terms 
in the perturbative series to all orders (exponentiation) 
that the vanishing of the cross section is recovered. 



2 Asymptotic Scaling 

Let us first sec what can be deduced on the fragmentation 
function from general considerations. If the quark mass m 
is infinitely large, hadronization effects are negligible, and 
the fragmentation function is just 5(1 — z). Taking a large 
but finite ratio m/A, one would expect the function to 
be somewhat smeared towards smaller z. This smearing is 
proportional to m/A, as expressed by the following scaling 
law (see e.g. 0): D(z) = (m/A)f((l - z)m/A). 

This property can be formulated more precisely upon 
taking moments, 



D(N,m 2 ) = I dz z jy ~ x D(z,m J ) 



N-l i 



(4) 



and it can be explicitly derived from the field-theoretic 
definition Q. One can consider two limits, one where the 
mass becomes large and the other where the moment in- 
dex N gets large. For large m one can match the matrix 
element © onto the heavy-quark effective theory, get- 
ting HI]: 



F(py,m 2 ) 

py 



exp (ipy) — ► T(j>yA/m) + (D(A/m), (5) 



namely, at the leading order in the large-m expansion the 
dependence on m and on the light-cone separation ?/_ (i.e. 
on py) is coupled: the matrix element becomes a func- 
tion of a single argument py A/m. Here A is the difference 
between the heavy-meson mass M and the heavy-quark 
mass m. For large N it follows from the definition J5J and 
from Q that 



D(N, m 2 



F(py, m 2 ) 



py 



exp (ipy) 



iN 



°^)- {(>] 



namely that to leading order in 1/N the JV-th moment 
of the fragmentation function can be obtained by analyt- 
ically continuing the matrix element as a function of the 
light-cone separation to the complex plane and evaluat- 
ing it at py — —iN. From and © together it follows 
that upon taking the simultaneous limit m — ► oo and 
N — > oo with a fixed ratio m/N, 



D(N, m 2 ) ~ T(pyA/m) 



iN 



so the fragmentation function becomes a function of a sin- 
gle argument NA/m. 

In Sec. 0] we shall see how the dependence on m and 
N through the combination NA/m follows from the large- 
order behaviour of the perturbative expansion in the large- 
st) limit. Having established Eq. Q non-perturbatively, 
we know that this is indeed the leading behaviour at large 
N and that corrections to this behaviour are suppressed 
by a power of 1/N. 

We see that the scale which characterizes the fragmen- 
tation process in the large z region is m(l — z) or, in mo- 
ment space, m/N. This scale has a clear meaning when 
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considering the bremsstrahlung off a heavy quark. Let us 
examine the emission in a frame where the quark energy 
E is much larger than its mass. The radiation pattern (to 
0(a s )) is 

dD C F q s 1 sin 2 9 

dzd sin 2 9 ~ 7T 1 - z J^T+T^/E 2 ) 2 ' 

where only the leading term in the limit z — > 1 was 
kept and the angle of emission 9 is related to the gluon 
transverse momentum by sin 2 9 = k\/ (E 2 z 2 {l — z) 2 ). As 
discussed in JOjj the radiation vanishes in the exact for- 
ward direction, but it peaks close to the forward direction 
at 9 ~ m/E (the 'dead cone'), or in a boost-invariant 
formulation at \k±\ ~ m(l — z). So m(l-z) is the typical 
transverse momentum of radiated gluons. The scaling law 
@ can be understood in physical terms as the observation 
that the hadronization effects (D(N,m 2 ) at large TV and 
m) are dominated by interaction with gluons of transverse 
momentum A — M — m. 



3 Factorization 

Factorization is based on the fact that dynamical pro- 
cesses taking place on well-separated physical scales are 
quantum-mechanically incoherent. This allows one to treat 
different subprocess independently of one another and to 
resum large corrections. 

Eq. © is often regarded as the separation between per- 
turbative and non-perturbative contributions to the cross 
section. However, factorization can be a much stronger 
tool upon considering separately the dynamics taking place 
on different physical scales. Consider, for example, the 
case of bottom production in e + e~ annihilation, shown in 
Figd Referring to @ one can naively interpret the gap 
between the data and some perturbative calculation as 
the "non-perturbative fragmentation function" and then 
try to bridge this gap using a model. As stressed above this 
interpretation leads to much confusion. Instead, the rea- 
sons for having large (perturbative and non-perturbative) 
corrections need to be identified and the corrections be 
resummed. 

The first step is to separate the scales involved. Upon 
neglecting higher order corrections which are suppressed 
by powers of m 2 /q 2 the moments of the cross section can 
be written as |11II12| 

&(N, q 2 ,m 2 ) = C(N, q 2 ;^ 2 F )E(N, ^ 2 F ,^ 2 F )D(N, m 2 ;f i 2 0F ). 

Choosing fi F ~ q 2 and /j,q F ~ m 2 , the coefficient function 
C and the fragmentation function D depend only on scales 
of order q 2 and m 2 , respectively. The evolution factor E 
can be obtained solving the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) equation. This factor then re- 
sums corrections depending on a s h\m 2 / q 2 to all orders. 
Resummation of this kind was implemented in computing 
the full line in Fig.^ Clearly, this is insufficient. 

Next, one observes that the subprocesses C and D may 
contain additional large corrections. One generic source of 



Jcl 



Y*(q) 




Fragmentation function 

Fig. 2. Factorization of bottom production in e + e" annihila- 
tion at large N. 

large corrections (see \1'S\) are running coupling (renor- 
malon) effects, which induce factorial growth of the coef- 
ficient at high orders owing to the increasing sensitivity 
to extreme ultraviolet or infrared scales. Infrared renor- 
malons in particular are non-summable and introduce a 
power-suppressed ambiguity in the perturbative defini- 
tion of any quantity. Since for observable quantities this 
ambiguity must cancel it can serve as a probe of non- 
perturbative contributions. 

Another source of large corrections develops at large 
N: the Sudakov logs |14U12j . As stressed above the frag- 
mentation process is dominated at large N by momenta 
of order m/N. When m and m/N become far apart the 
concept of factorization applies again, and can be used to 
resum logs of N into a Sudakov form factor. This resum- 
mation takes the form of exponentiation in moment space. 
A similar situation occurs in the coefficient function C, as 
is demonstrated in Fig. [3 C is dominated at large N by 
the invariant mass q 2 /N of the unresolved jet which recoils 
against the measured heavy meson. The fact that this jet 
was also initiated by a heavy quark plays no role at this 
level : the relevant scale here is the total invariant mass 
of the jet. The same jet function dominates deep inelastic 
structure functions at large N [15116) . 

It should be emphasized that factorization (contrary 
to its diagrammatic proofs) is a non-perturbative concept. 
One should therefore expect that non-perturbative correc- 
tions on a certain scale would factorise together with the 
corresponding perturbative sum. In particular, this must 
apply to renormalon-related power corrections. In the case 
of Sudakov logs factorization leads to exponentiation. Go- 
ing beyond the logarithmic level, one finds that power cor- 
rections on the corresponding scale exponentiate as well. 
This is the conceptual basis for the "shape function" ap- 
proach to hadronization corrections, which has been de- 
veloped in the context of event-shape distributions [17U18I 
IT§] (see also [20] )• This is also the basis of the approach 
of |15Ulfi| to higher twist in deep inelastic structure func- 
tions at large N and of our approach to heavy-quark 
fragmentation. 



4 Dressed Gluon Exponentiation 

In order to deal with heavy-quark fragmentation at large N 
both Sudakov logs and renormalons need to be taken into 
account. At large N, the perturbative coefficients are dom- 
inated by Sudakov logs. However, the resummation of the 
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P/z 



Off-shell Quark 




On-shell Quark 

Fig. 3. The diagram contributing to the fragmentation func- 
tion to leading order in the flavour expansion (the large-/3o 
limit) in the A ■ y = gauge. 

leading logarithms alone does not provide any information 
on power corrections. It is the subleading logs generated 
by the running of the coupling which produce the renor- 
malon ambiguity |19U21llT] . Their resummation is there- 
fore essential to probe the non-perturbative regime. 

From these considerations it follows that the Sudakov 
exponent needs to be computed to all orders rather than to 
some fixed logarithmic accuracy. Clearly, the full calcula- 
tion cannot be done. However, relevant all-order informa- 
tion can be obtained from the large-/3o limit correspond- 
ing to a single dressed gluon. Calculating the Sudakov 
exponent in this way is referred to as "Dressed Gluon Ex- 
ponentiation" (dge) poHirr] . 

A process-independent calculation of the fragmenta- 
tion function Q in the large-/?o limit was performed in pQ. 
In the light-cone axial gauge A ■ y = where the path- 
ordered exponential is 1, there is just one diagram - see 
Fig-El This diagram was computed using an off-shell gluon 
splitting function, which was derived identifying the limit 
where the massive quark propagator prior to the emission 
of the gluon is singular 1 . 

The result for the logarithmic derivative of the frag- 
mentation function, written as a scheme invariant Borel 
transform, is: 



dD{N,m 2 ) C F [°° J [A 2 

z 



Q3' 



dz i 



.JV-1 



(1-zf 



1 



I z (l-u) + -(l-z)(l + u) 



(9) 



where A is in the MS scheme. A generalization of this 
result beyond the large (3q limit which fully captures the 
next-to-leading logarithms (NLL) was constructed in pQ . 

Eq. takes into account the cancellation between real 
(z N ~ v ) and virtual (1) corrections. In the square brackets 
we distinguish between z = 1 singular and regular terms. 
The former lead to logarithmically enhanced contributions 



1 In this limit the gluon virtuality k 2 
turn k 

neously keeping the ratios between them fixed. This is a gen- 
eralization of the quasi-collinear limit discussed in |22I12| . 



its transverse momen- 
2 and the quark mass m 2 are taken to be small simulta- 
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Fig. 4. The perturbative DGE result compared with ALEPH 
data for B production in e + e _ annihilation, plotted as a func- 
tion of the moment N. 



solid: DGE with Principal Value 
dashes: NLL 

dots: no soft-gluon resummation 



pJii 



T 



TT 




Fig. 5. Same as Fig. 2] but plotted as a function of the energy 
fraction of the measured heavy meson. 

in the perturbative expansion, and therefore need to be 
exponentiated. 

According to © the natural scale for the renormal- 
ization of the coupling at fixed z is (1 — z) 2 m 2 / z. Thus, 
integrating over the Borel variable u first is not possi- 
ble for (1 — z)m < A. As expected, perturbation theory 
breaks down when the gluon virtuality or its transverse 
momentum become comparable to the QCD scale. This 
constraint takes a completely different form when consid- 
ered in moment space: infrared renormalons show up. 

We proceed to compute the Sudakov exponent in the 
large-/3o limit by isolating the z = 1 singular terms, per- 
forming the z-integration and then integrating over m 2 . 
The result is: 



, 2 2 s Cf f°°du [A 2 



(10) 
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2 \ u 

to \ 



B A {u) \nN - Bl GE (u,N) 



where 



nDGE 

B D 



(u,N) 



;* u (1 - u) r(-2u) (N 2u - 1) . (11) 



The to 2 integration requires to introduce an ultraviolet 
subtraction: a /i^ -dependent counter term which cancels 
the u = singularity of the fragmentation function. This 
term is the well-known cusp anomalous dimension [231124) . 
given by B^(u)lnN, where B^(u) = 1 + |u + . . . (we 
use the MS factorization scheme). Note that contrary to 
B^ E (u, N) this subtraction term has just a single In N to 
any order in u and it is also free of infrared renormalon 
singularities. 

According to Eq. (fTT|) , renormalons in the Sudakov ex- 
ponent (|10|) appear at all integer and half integer u values 
with the exception of u = 1. It is clear from Eq. JHJ that 
these renormalons are exclusively related to the z — > 1 
limit. To define the perturbative sum corresponding to 
\nD(N, m 2 ) one needs to integrate over u with some pre- 
scription that avoids the poles. The natural choice is the 
principal-value (PV) prescription (it was implemented nu- 
merically in pP). The ambiguity in choosing a prescription 
is compensated by power corrections corresponding to the 
residues. Introducing a free parameter for each singularity 
one ends up with an additive correction to the perturba- 
tive Sudakov exponent having the form: 

NA fNAV [NA\ a 

lnD NP (NA/m) = -e x e 3 -£4 +■ • • ■ 

to \ m J \ m J 

(12) 

Finally exponentiating the result to compute D{N,m 2 ) 
the perturbative and non-perturbative contributions ap- 
pear as two factors: 

D(N, m 2 ;^ 2 F0 ) = D PT (N, to 2 ; /j| ) D NP {NA/m). (13) 

The leading power correction of the form NA/m predicted 
in is readily obtained from (|12f> upon expanding the 
exponent. 

It should be stressed that in both the perturbative 11011 
and the non-perturbative (|12[> contributions to the frag- 
mentation function we considered here only the leading 
terms at large N. At the perturbative level the result is 
improved pQ by matching it with the full NLO coefficient. 
At the non-perturbative level, there may be additional 
0(A/m) terms which we do not parametrize, and con- 
sequently the description of the first few moments is of 
limited accuracy. In practice, to deal with low moments, 
it is useful to modify the parametrization l|12)l replacing 
N — > N — 1 such that the N = 1 moment is exactly 1, 
as it must be by definition. 

The perturbative PV-regulated DGE result of Eq. ffT77f> . 
matched to the NLL and the NLO expressions and com- 
bined with the proper coefficient function, is compared 
as a function of N with the ALEPH data in Fig. ||| (full 
line) . In contrast with the NLL result of Ref . ^21 (dashed 
line), the DGE one does not have a Landau singularity 





1 " " 1 " 1 'E 


\ e, = 1.23 


± 0.23 


- \ A = 0.217 


± 0.025 GeV - 


fit range ^Jjj 




; , , , 1 , , , , 1 , , , 





10 20 30 40 50 

N 




Fig. 6. A DGE-based fit to ALEPH data. The free parameters 
are a a (i.e. A) and the leading power correction coefficient ei, 
which controls the shift of the perturbative distribution. 

and thus it extrapolates smoothly towards the values of 
N <; to/ 'A which are beyond perturbative reach. 

Also shown in Fig. 01 is the ambiguity (band shown by 
two dot-dashed lines) corresponding to the residue of the 
first renormalon pole located at u = 1/2. The lower edge of 
the band just matches the data, indicating that the power 
correction of the form and magnitude (!) expected based 
on the renormalon analysis is supported by the data. 

The different perturbative results are converted to x 
space in Fig. [3] Here the significant impact of Sudakov 
resummation to NLL as well as that of the additional 
renormalon resummation achieved by DGE on the shape 
of the distribution is evident. Note that the shape of the 
DGE curve resembles that of the data but it is centered at 
larger x. Indeed, the leading effect of the non-perturbative 
function (assuming in l|12(l that only t\ ^ 0) is a shift of 
the entire perturbative distribution, very much the same 
as the leading corrections in event-shape distributions [171 
I20II19| . Finally, regarding the non-perturbative parame- 
ters e n as free parameters in a fit, the data can be well 
described. The result of a fit in moment space where the 
only non-perturbative correction is e% is shown in Fig. El 

Upon using more non-perturbative parameters the de- 
tails of the prediction i|12|) can be confronted with data. 
The analysis in £Q shows that subleading non-perturbative 
corrections at the exponent are rather small, and the ab- 
sence of a correction of the form N 2 A 2 /m 2 can be consis- 
tent with the data. 



5 Conclusions 

We described here a new approach to the QCD descrip- 
tion of heavy-quark fragmentation concentrating on the 
z — ► 1 limit. It was first rigorously demonstrated that 
the non-perturbative dynamics is dominated by the scale 
to(1 — z). This scale corresponds in perturbation theory 
to the transverse momentum of gluons radiated from the 
heavy quark. Based on a renormalon analysis we extended 
the perturbative technique for resumming soft gluon radi- 
ation to the non-perturbative regime, identified power-like 
effects and separated them from the perturbative fragmen- 
tation function by means of a PV prescription. The non- 
perturbative contribution was then parametrized based on 
the renormalon ambiguity. We found that the simplest 
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possible parametrization of power corrections which fol- 
lows from renormalons, namely a shift of the perturbative 
distribution, is sufficient to describe the data on B pro- 
duction in e + e~ annihilation. This way phenomenological 
models for the non-perturbative fragmentation function 
are not needed. 

The fragmentation function was treated, based on its 
definition (TjQ), in a process independent way. The results 
are thus applicable independently of the production pro- 
cess, given that the corresponding coefficient function in 
the MS scheme is known. Universality of the leading power 
corrections at large z can now be tested experimentally. 
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